The Condition for Universality at Resonance and Direct Measurement of Pair 

Wavefunctions Using rf Spectroscopy 
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We show that when the Fermi energy of a Fermi gas is much smaller than the intrinsic energy width 
of a Feshbach resonance, the system behaves like a Fermi gas with a contact potential. This in turn 
implies universality at resonance, and large fermionic pairs in the strongly interacting regime. Recent 
experiments of JILA 1] and MIT 2] turn out to be deep inside the universal regime, which explains 
the perfect fit of these experiments by the BEC-BCS crossover theory with contact potential^. We 
also show that rf spectrocopy can be used to map out the pair wavefunction directly. 
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It is by now well known that a weakly interacting 
quantum gas can be made strongly interacting by tun- 
ing the system close to a scattering resonance, where the 
s-wave scattering length becomes divergent. There are 
many ways to generate a scattering resonance. The sim- 
plest one is to vary the pair potential on the microscopic 
scale. This, however, is difficult to achieve. An alternate 
way, which leads to the recent explosion of activities, is 
to use Feshbach resonance. This resonance is achieved 
by Zeeman shifting the energy of a bound state in a 
closed channel to zero energy, thereby generating con- 
siderable resonance scattering for particles in the open 
channel. In the last fourteen months, many remarkable 
properties of Fermi gases near a Feshbach resonance were 
discovered - universal interaction energy 0, Q, molecu- 
lar condensates p, and the long sought condensation of 
fermion pairs. Evidence of such condensation was re- 
ported by C. Regal, M. Greiner, and D. S. Jin 1] three 
months ago in 40 K, and a month ago by Ketterle's group 
in 6 Li[2|. Evidence of superfluidity has also been ob- 
served recently by Grimm's Q and Thomas' Q groups. 

As these phenomena emerge, there are questions about 
whether they are general properties of scattering res- 
onances, or specific properties of Feshbach resonances. 
This question is sometimes phrased in a narrower con- 
text as whether single channel models and two-channel 
(or resonance) models have the same physics near reso- 
nance. The former refers to models with two types of 
fermions interacting with a pair potential. The latter 
are those that describe open channel and closed channels 
physics explicitly. Another important question is under 
what conditions are the properties of these models uni- 
versal, i.e. independent of microscopic details. Moreover, 
if both types of models are universal, do they belong to 
the same universality class? These questions are not only 
of theoretical interest, but also of practical importance in 
interpreting current data and guiding future experiments. 

The purpose of this paper is to point out the condi- 
tion for universal behavior (referred to as "universality 
condition" ) for both single channel and two channel sys- 
tems, its implications on the ground states, and a simple 
method to measure the pair wave function directly. We 
shall show that (1) If the Fermi energy is well within 



an intrinsic width of the resonance (defined later), then 
the interaction becomes <5-function like, which guaran- 
tees universality at resonance. (2) There is only one 
universality class, given by the behavior of Fermi gases 
with contact interaction. (3) Universality implies that 
the condensate must have large pairs (of the size of inter- 
particle spacing), where closed channel bound states play 
an insignificant role. (4) The recent JILA and MIT ex- 
periments are deep in the universal regime. This justi- 
fies the use of a single channel approach for the JILA 
and MIT experiment and explains the perfect agree- 
ment between the phase boundary in ref . 0, Q and the 
T c predicted by the crossover theory with 5-function 
potential3]. (5) Condensates with molecular rich (and 
hence small) pairs can be found in "narrow" resonances. 
(6) RF spectroscopy can be used to map out the pair 
wavefunction across the resonance directly. 

The emergence of universality is most obvious for 5- 
function potentials (in single channel models) with cou- 
pling constant g = Aith 2 a s /M , where a s is the s-wave 
scattering length and M is the mass of the fermion. At 
resonance, a s diverges. The only remaining length scale 
is then n" 1 / 3 , where n is the density. As a consequence, 
the thermodynamics of the system becomes universal 0, 
i.e. independent of microscopic details. It is clear that 
universal behaviors only emerge when all other length 
scales become irrelevant. However, all resonances have 
an intrinsic length scale r* , which is related to the width 
of the resonance in energy space. To be precise, con- 
sider the scattering length near a Feshbach resonance, 
a s = a,b g (l — W/(B — B )), where a bg is the background 
scattering length, W is the width of the resonance, and 
B Q is the location of the resonance. Near resonance, it is 
sufficient to focus on the resonance term and write 

a s = -ti/v, v = Hco(B - B ), r\ = a bg fi co W (1) 

where /i co is the difference in magnetic moment between 
fermions in closed and open channel, and v is the de- 
tuning energy. On can then define an intrinsic length 
scale r* or intrinsic width (A_B); n in magnetic field as 
h 2 /2Mr* =n = r*fi co (AB) in , or 
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In turn, one can write 

a s /r* = (A-B)in/ (B Q - B) (3) 

[For single channel systems where a bound state is gen- 
erated when an energy parameters u (such as the depth 
of the well) reaches a special value u a , a s can still be 
written in the form eq.Q with v = u Q — u and 77 be- 
ing some constant]. The presence of r* means the many 
body system has an additional dimensionless parameter 
y = kpr*, where kp = (S^tt,) 1 / 3 is the Fermi wavevec- 
tor. We shall show that reduction to contact potential 
(hence universality at resonance) occurs when 

y = k F r* « 1, or fi co (AB) in / Ep » 1. (4) 

Other related works: Before proceeding, we would 
like to mention that eq.© has also been recently identi- 
fied by Eric Cornell as the condition for the equivalence 
between two channel and single channel models, (see later 
discussions) . While we agree with his conclusions, we 
point out eq.© is in fact a (stronger) condition for uni- 
versality, which implies channel equivalence. Such dis- 
tinction is important because not all single channel mod- 
els are universal, as we shall see later. (The different po- 
tentials chosen in ref.Q happens to be all in the universal 
regime). In a recent prepint, G. Brunnflll] has also iden- 
tified eq. © as the condition for universal behavior when 
considering a normal state at T — without superfluid 
correlations. The lack of superfluid correlation makes it 
difficult to connect to current experiments. During the 
writing of this paper, a preprint by De Palo et al.[l2^ 
has appeared in which they come to similar conclusions 
for the single channel problem. The present work differs 
from those mentioned above not only in perspective and 
approach, but also in the conceptual emphasis of under- 
lying importance of universality, which controls channel 
equivalence and the nature of the pairing state. In ad- 
dition, we point out how rf spectroscopy can be used to 
probe the pair wavefunction in the universal regime. 

Before ending this section, we would like to para- 
phrase Cornell's arguments[10|, which we find illuminat- 
ing. They should be appreciated together with the re- 
sults of the explicit calculation presented later. Cornell 
notes that (i) For a two body system, when \B — B Q \ << 
(A£?)i n , or a s >> r* in eq.©, the population of the 
closed channel is very small, and the problem should 
therefore be single channel like . (ii) For a Fermi gas 
with density n, as B approaches B Dl two things can hap- 
pen. Either one first enters the "single channel region" 
\B - B D \ « (AB) (hence a s /r* >> 1) and then the 
strongly interacting regime n|a s | 3 >> 1, or vice versa. 
In the former case, the system become single channel 
like before becoming strongly interacting, (iii) The con- 
dition for being single channel like can be obtained by 
replacing v in eq.© by the Fermi energy Ep, which is 
the "detuning" of the fermions at the Fermi surface. The 



"single-channel" condition is then ii co {AB) la / 'Ep » 1, 
which is eq.©. We shall now prove the condition eq.© 
using the two-channel model. 

Universality condition for two-channel mod- 
els: We shall consider the resonance model 14] H — 
= Ek>k - A*)«L«k,. +£k(ek/2 ~2[i + V)blb k 

+ rl/2 «Ek ! ,( i !, 8 k +q /2/-k+,/2j + h - c -) where a L 

creates a fermion with momentum k and spin a, {a =] 
,|); b\. creates a closed channel bound state with mo- 
mentum k, £1 is the volume of the system, and a is the 
(intensive) coupling between the closed channel bound 
state and the open channel scattering state. We shall 
ignore the background scattering length term since it is 
unessential for the resonance physics. The quantity V is 
the "bare" detuning which is related to the physical de- 
tuning by an infinite constant T> = v + ^- ^ k so as 
to cancel any unphysical ultra-violet divergences in the 
problem. With this normalization, it is straightforward 
to solve for the two-body T-matrix and find the s-wave 
scattering length, which is a s = — ^7- 

Next, we study the ground state properties of the res- 
onance model using the BEC-BCS crossover theory [lij 
where one assumes a condensation in the bound state 
(& q=0 ) = $ m in the closed channel and a condensation 
of zero momentum pairs (ak,fQi-k,j.) m the open channel. 
Standard approach then gives the mean field equation 

= ^{2^^ 2^)' (5) 

where Ek — yj (et A*) 2 + A 2 , and A is the gap param- 
eter related to $ m as A = a^/n mi n m = |<i> m | 2 /f2, where 
n m is the density of closed channel molecules. The chem- 
ical potential is also constraint by the number equation 
n = 2n m + i Ek, CT ( a k, CT a k, CT >> or 

k v 7 

Eq. (J5J and © together give fi and A (hence the molecu- 
lar fraction n m /n in the closed channel) as a function of 
n and the physical detuning v for a given system (a). 

It is useful to write eq.© and © in dimensionless 
form. Defining q = k/kp, jl — /i/Ef, A = A/Ep, where 
Ep = h 2 k F /2M; and noting that n = kp/Sn 2 , it is easy 
to show that eq. © and © can be written in dimension- 
less form y{u-2p) =Fi(p,,i>), (4/3tt 2 ) = yA 2 +F 2 (/i, A), 
where F\ and F2 are functions of jl and A|l||. Micro- 
scopic details are characterized by y only. 

The emergence of universality is best illustrated in fig- 
ure 1 where we have plotted jl at resonance v = as a 
function of 1/ ^fy . (There is no particular significant in 
this choice of plotting variable except it gives a clear dis- 
play of approach to universality). We see that as y^ 1 ^ 2 
exceeds beyond 10, Jl becomes independent of y (hence 
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FIG. 1: Chemical potential on resonance (y = 0) as a function 
of y- 1 ' 2 , y = kpv* . The solid line and circles are results 
of two-channel and single channel calculations respectively. 
Both models exhibit universality (independent of microscopic 
detail y) for i/ -1 ^ 2 > 10. The corresponding values in the 
JILAQ and MIT@ experiments are y~ 1/2 = 15 and 80, which 
are inside and deep inside the universal regime. 

microscopic details). Moreover, ft saturates at 0.59, the 
same value of the single channel <5-function potential. At 
the same time, the fraction of close channel molecules 
n m /n = \$ m \ 2 /N, (n m = A 2 /a 2 ) drops rapidly as y~ 1/2 
increases, as shown in figure 2. This shows that the molec- 
ular component is insignificant in the universal regime. 
The size of the fermion pair must therefore be of order 
n -i/3_ rjv^g k enav i or f ft across resonance (i.e. plotted 
for different l/(kpa s )) for different a is shown in figure 
3. For y~ x l 2 > 10, all ft coincide, and is precisely that 
given by single channel systems with (^-function potential, 
demonstrating the universality of the contact potential in 
the wide resonance regime. 

The reduction to single channel ^-function can in fact 
be proven simply as follows. By noting that the term 
v/a 2 in eq.© is simply —1/g, where g = 4irh 2 a s /M , 
eqs.© and © reduce to 

~~g = 2^)' (7) 
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FIG. 2: The fraction of closed channel molecules on reso- 
nance v = as a function of j/ _1//2 , y = kpr* . One sees that 
the fraction of closed channel molecules is insignificant in the 
universal regime. 

values of y = kpr* are generated. In figure 1, we have 
also displayed (with circles) the chemical potential along 
a path in parameter space which gives a monotonically 
increasing (There are many such paths). One 

sees from fig.l that the chemical potential again reaches a 
constant when y _1//2 >10, which is a sign of universality. 

Relation to current experiments: In the JILA 
experiment Q], one has n = 10 13 cm~ 3 , W = 8 G, 
ab g = 170a b, t^co ~ 2^s, where as and [Ib are the Bohr 
radius and Bohr magneton. Since r\ = ^ C oWab g , we find 
from eq.© that y~ 1/2 = 15, (or y = k F r* = 0.004) 
which is inside the universal region as shown in figure 
1. In the MIT experiment 0, one has kp 1 = 2000as, 
W = 180 G, ab g = — 2000a_B, fi c0 ~ 2/is, which implies 
j/ -1 / 2 ~ 80, (or y = 1.5 x 10~ 4 ), which is deep inside 
the universal single channel regime. This also explains 
the surprisingly good agreement between the prediction 
of crossover theory using a (5-function potential with the 
experiments in ref . and 0] as pointed out in ref. Q . 

Direct Measurement of Pair Wavefunction: We 
conclude by pointing out that rf spectroscopy can be 



(8) 



in the limit of large a 2 since n/a 2 and A 2 /a 2 can be 
dropped from eqs.© and © respectively. Eq.Q and 
(JSJ) are precisely the gap equation and number equation 
for the single channel systems with contact potential. 

Single channel models: It is important to note that 
not all single channel models are universal. This is also 
pointed out in ref. 12]. To show that the condition for 
universality is also given by eq.Q), we consider a square 
well potential with a barrier: V(r) = —u < 0, r < r ; 
V(r) = mi > for r„ > r > n, and V(r) = for 
r > r\. As the parameters u Q and u\ are varied, different 
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FIG. 3: The behavior of chemical potential across resonance 
for different values of y _1,/2 : 0.25 (dot-dashed), 1 (dotted), 
5 (solid grey), 10 (dashed), 25 (solid black). For y > 10, all 
curves coincide with that of single channel contact potential. 
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and reconstruct the real space wavefunction. The behav- 
ior of ea. l|10|) for fermions interacting with a <5-function 
potential in BCS, resonance, and BEC regimes are shown 
in figure 4 for the case of W > 0. The case of W < 
can be obtained by flipping about the vertical axis, i.e. 
changing hui — W to W — tko. Note that the shape of 
the curve in the molecular regime is consistent with the 
shape observed in ref.0] (which has W < 0). 

We thank Hans Peter Buchler for a discussions on rf 
spectroscopy. This work is supported by NASA GRANT- 
NAG8-1765 and NSF Grant DMR-0109255. 



FIG. 4: Transition rate in rf spectroscopy experiments (in 
arbitrary units) as a function of frequency. The values of 
(fcfO,,)" 1 for different curves are: 2 (dashed grey), 1 (dashed 
black), (solid black), -1 (solid grey), -2 (dotted). 



used to determine the density profile of the pair wave- 
function directly. It will therefore be a useful tool to 
explore the properties in the universal regime |l8j|. Con- 
sider exciting a fermion to a different atomic state 
c. For example, aj and c can be the hyperfine states 
F z = —5/2 and —7/2 respectively. Let W be the en- 
ergy difference between aj and c. The inclusion of the 
c particle (of the same mass) and the rf field will in- 
troduce an additional term ^ k £fcC k c k to the unper- 
turbed Hamiltonian and a time dependent perturbation 
V = A^ k (c k a k ^e~ luJt + h.c). The rate of transition 
from ground state to the c state is given by the Fermi 



Golden Rule, TZ = 2zl (/|Ac k a k)T |G) 5{£ f -£ t -huj), 

where \G) and |/) are the ground state and the excited 
states of the system, with energies £, t and £ f respectively. 
Considering an initial fermion ground state with 2N par- 
ticles with energy E^ N , we have £{ — E^ N . Since the ex- 
cited state consists of a Bogoliubov particle with energy 
i?k = \J (ek — n) 2 + A k on top of a ground state with 
energy £^v-i an d a particle c k with energy ek + W, we 
haveft(w) = ^iEk^C^k + ^ + ek-M-M. OT 



K(u) 



2ir\ 2 1 



D(e*) 



1 + dE k /de k 



(9) 



where = \/[l ~~ ( £ k n)/Eu)]/2 is the coherence 
factor, D(e) — (27r) _3 <i 3 fc/dek is the density of state, 
and e* is the solution of the equation hco — W — 
y/ (e* — /i) 2 + A 2 + e* — [i. For (5-function potentials, Ak 
is a constant, we then have I + dE^/de^ = 1 — u k = it k . 



Eq.© then reduces to 



(10) 



Since v^/u^ is the Fourier transform of the pair wavefunc- 
tion apart from a normalization constant [15|, and since 
D{e) on a/c; one can therefore extract directly from the 
signal of rf spectroscopy fea. Hl(J|) ) the coefficient i>k/wk 
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